Two scenarios for avalanche dynamics in inclined granular layers 
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We report experimental measurements of avalanche behavior of thin granular layers on an inclined 
plane for low volume flow rate. The dynamical properties of avalanches were quantitatively and 
qualitatively different for smooth glass beads compared to irregular granular materials such as sand. 
Two scenarios for granular avalanches on an incline are identified and a theoretical explanation for 
these different scenarios is developed based on a depth-averaged approach that takes into account 
the differing rheologies of the granular materials. 
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Avalanche behavior of granular material has many nat- 
ural realizations from snow avalanches to massive rock- 
slides. Physics-oriented avalanche investigations focusing 
on sand-pile kinematics [l| have given way to more recent 
work, which has probed granular interactions and the 
appropriate balance between continuum and discrete ap- 
proaches |2|. Often avalanche dynamics have been stud- 
ied with a bulk granular substrate, either in a rotating 
drum [D, 0, 0, IS wnere the critical angle C , i.e., the 
angle where grains start to flow, is reached by the slow ro- 
tation of the drum, or on a pile [3, 0] • Alternatively, one 
can investigate flow on an inclined plane where an under- 
lying solid surface constrains the flow. In inclined layer 
granular flow at high volume flow rates a uniformly thick 
layer forms, whereas for lower flow rates waves in the 
form of thickness variations appear 0, 0, 0, 0, 0, ^( . 
For both the freely flowing and wave modulated cases, 
continuum descriptions of the flows based on flow rheol- 
ogy describe the experiments well 0, 0, ^3 • 

In this Letter we focus not on steady or modulated 
flows but rather on flows exhibiting distinct, well-resolved 
avalanches. We have explored these flows on an inclined 
plane and have discovered two distinct scenarios for the 
dynamic avalanche behavior, depending on the character 
of the grains. For rough non-spherical grains (RNSG), 
avalanches are faster, bigger and overturning. Individ- 
ual grains have down-slope speeds that exceed the front 
speed. By contrast, avalanches of spherical glass beads 
(SGB) are quantitatively slower and smaller, and the par- 
ticles always travel slower than the front speed. We show 
that this difference in behavior arises from the differing 
rheologies of the different particle types in steady-state 
flow, and is linked to the stability of such flows. A theory 
based on the nonlinear properties of the depth-averaged 
equations makes quantitative predictions that agree with 
the experimental results and suggests that our different 
avalanche structures are related to the difference between 
roll waves and flood waves in hydrodynamic contexts [HI . 
In particular, the propagation of the RNSG avalanches 
involves a shock, while the SGB avalanches have a con- 



tinuous structure. 

Two key lessons emerge from this study. The first is 
that our understanding of the rheology of dense gran- 
ular flows is now sufficiently robust to allow successful 
modeling of quite detailed dynamical phenomena such 
as avalanche profiles. The second is that the quantita- 
tive differences between the rheology of the most com- 
monly studied spherical particles and the more realistic 
rough particles result in qualitative, indeed, dramatic dif- 
ferences in their dynamical behavior. 

A thin layer of granular material inclined at an angle 6 
is stable for a wide range of and thickness h. The onset 
of flow is expected only above a critical layer thickness 
[20L EH /i c , and the flow subsides at h s < h c . The val- 
ues of h s (0) and h c (0) decrease rapidly with increasing 
6. Thus, the layer can become unstable by increasing 
by a small amount SO |J [H| or by increasing h by 
adding new grains to the layer at a low flow rate, as 
we do here. In the former case, where the whole layer 
becomes metastable 0, [HI , the shape and propagation 
of the avalanches depends critically on 66 such that ei- 
ther strictly downward or simultaneous downward and 
upward expanding avalanches can occur depending on 
the magnitude of SO. In our case a homogeneous static 
layer is prepared that is stable to small perturbations. 
New grains are only added at the top region (5%) of the 
plane in the form of a very low incoming flux (shower) 
perpendicular to the plane. Due to the low incoming flux 
the height slowly increases locally, and when the little pile 
formed in this manner becomes unstable an avalanche is 
created, which travels down the rest of the plane on top 
of the stable static layer. 

The experiments reported here used an inclined plane 
2.2 m long and 0.4 m wide with a rough surface of parti- 
cles glued to a glass substrate or sandpaper with different 
grit sizes. The grains were of two types: RNSG of sand 
or salt with several mean diameters d (sand: 0.4±0.1 
mm; sand: 0.2±0.1 mm; salt: 0.4±0.1 mm) and SGB 
with d = 0.5±0.1 mm. The mass flow rate per unit 
width of the channel Q was adjusted so that individual 
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well-resolved avalanches formed, and was 0.17 g/s-cm for 
RNSG and 0.05 g/s-cm for SGB. Varying the incoming 
flux or the kinetic energy of the incoming particles gives 
rise to to a change in the size distribution and frequency 
of avalanches, which is not the subject of the present 
study. The angle 9 was varied between 32° and 41° for 
RNSG and between 22° and 26° for SGB. The differences 
in Q and for the different materials reflect the material 
variations in the critical angle 6 C and the angle of repose 
r where flow stops. The experimental methods for ob- 
taining the data presented below include: visualization 
of grain motion with high speed (lOOOfps) video imag- 
ing, determination of the lateral sizes of avalanches by 
image differencing, reconstruction of the 2D height pro- 
files of avalanches using a laser sheet, and determination 
of particle and front velocities by the analysis of single 
particle trajectories on space-time plots. 




FIG. 1: Height profiles of a) RNSG avalanche for = 36.8°. 
Image size 7.2 cm x 56 cm, (vertical size rescaled by 25x) 
maximum height: h m = 0.34 cm, static layer thickness h s = 
0.12 cm; b) SGB avalanche for = 24.3°. Image size 12.2 cm 
x 46.8 cm, (vertical size rescaled by 25x) maximum height: 
hm = 0.29 cm, static layer thickness h s = 0.18 cm. 

Figure Q] shows the reconstructed 2D height profiles for 
sand and glass bead avalanches for operating conditions 
(6 = 36.8°, 24.3° and Q = 0.17 g/s-cm, 0.05g/s-cm for 
RNSG and SGB, respectively) that produce roughly the 
same number of spatially-localized avalanches in an in- 
stantaneous image. For gravity- driven flows on an in- 
cline 22] the characteristic length scale is the height 
of the layer h s and the corresponding velocity scale is 
y/gh s cos 6. This normalization collapses data for differ- 
ent sized RNSG as shown in Fig. [2| where the spread 
is about 25% around the mean of the three data sets. 
For the SGB the velocities of avalanches with the same 
dimensionless areas are smaller by roughly a factor of 
4. Further, the maximum dimensionless area of SGB 
avalanches is less than for RNSG avalanches, again by 
about a factor of 4. The corresponding dimensionless 
avalanche velocity as a function of avalanche height, nor- 
malized by h s , is shown in Fig. 03. The data support 
the hypothesis that there are two classes of avalanches: 



SGB avalanches have maximum height h m that is always 
less than 2h s whereas the height of RNSG avalanches is 
always greater than 2h s . Our measurement of the height 
of avalanches relative to h s for SGB h m /h s w 1.45 ± 0.1 
is very similar to previous measurements for granular 
waves |9| for which h m /h s « 1.55 ±0.1, whereas the 
RNSG heights reported here are considerably higher: 
hm/h 8 « 2.5 ±0.2. 
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FIG. 2: Dimensionless avalanche front velocity vs a) dimen- 
sionless area A and b) dimensionless maximum height hm/h s 
for sand (x), glass beads (o), salt (*) and fine sand (A). 

The origin of this dramatic difference is revealed in 
detailed measurements of the velocity distribution of 
surface particles for the two types of avalanches. The 
avalanche front speed Uf and the surface mean parti- 
cle speed u p right behind the front can be determined 
by measuring the image intensity along the centerline 
of an avalanche. Plotting the intensity in a space-time 
plot, see Fig. 03 reveals streaks associated with particles 
and a well-differentiated front for each avalanche type. 
For RNSG avalanches (Fig. the streaks within the 
avalanche are less steep than the front line, indicating 
Up > Uf, whereas the opposite is true (Fig. Eb) for the 
SGB avalanches. In Fig. 21 the plot of u p versus Uf shows 
that for SGB u p < Uf, with a ratio u p /uf « 0.7, whereas 
for RNSG avalanches u p > Uf with u p /uf « 1.4. An- 
other feature of the avalanches seen in Fig.[3|is the contin- 
uous form of SGB avalanches as evidenced by the curved 
paths just ahead of the front (representing the accelera- 
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tion of particles in this region) as compared with particles 
being thrown out of the main body of the avalanche for 
the RNSG type. High speed imaging demonstrates [2^ 
that an RNSG avalanche consists of a fast moving packet 
of grains that overtakes the front, like a breaking wave 
in a fluid, which entrains granular material from the un- 
derlying layer as it passes over. SGB avalanches, on the 
other hand, are continuous in that the static stresses that 
hold grains in place in front of the granular packet change 
and the previously stable packing collapses. 
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consider the rheology. For steady- state flows, both sand 
and glass beads have depth- averaged velocities u given 
by the Pouliquen form |2(J, |22|, |24| 



F(h/h a {0)), where F{y)=/3y-j, (1) 



gh 



where h s {6) is the minimum height of a flowing pile for 
an angle 0. (Note that the depth- averaged quantities 
of the continuum model are different than the surface 
velocities measured in the experiments.) For sand [Io( . 
the coefficients in the function F(y) are (3 « 7 « 0.7, 
whereas for glass beads, /3 ~ 0.14 and 7^0 0, H3, 
Eij . We now apply these results to our RNSG and SGB 
systems, respectively. 

The simplest dynamical picture of granular flows down 
an incline is achieved using the Saint- Venant shallow flow 
equations, adapted for granular media by Savage and 
Hutter [25|. For a flow of height h and mean velocity 
u, these describe flow down a plane, with the plane par- 
allel to the ^-direction, by the averaged equations 



dh d(hu) 
dt dx 



(2) 



FIG. 3: Space- time plots taken along the symmetry axis of 
the avalanches a) sand for = 36.8° with scaled velocities 
(Froude numbers) Fr/ = 2.12, Fr p = 2.55, and b) glass beads 
23.3° with Fr/ = 0.55, Fr p 0.33. 
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FIG. 4: Ratio of particle and front velocities u p /uf for glass 
beads (o) and sand (x). The inset shows u p vs Uf for = 36.8° 
(sand) and = 25.2° (glass beads). 

These two distinct types of avalanches can be under- 
stood as consequences of the underlying rheology of the 
RNSG and SGB systems and of the structure of the dy- 
namical equations for a propagating front. Let us first 



d(hu) d(hu 2 ) 



dt 



dx 



tan# — fi(u, h) 



K — — ghcost 
ox J 



(3) 

Here a is determined by the profile of the flow, a = 1 for 
plug flow (as in Ref. |25j). a = 4/3 for a linear flow pro- 
file, or 5/4 for a convex Bagnold profile The param- 
eter K is determined by the ratio of the normal stresses 
in the flow: the stress parallel to the bed, <j xx , and that 
perpendicular to the bed, a zz . Numerical results show 
that K = cr xx /(j zz ~ 1 for steady-state flows 24]. The 
friction coefficient h) is determined by the require- 
ment that the steady flow obey the rheology shown in 
Eq. fQ, and will thus vary with the particle type. 

Since the dimensionless velocity, typically given as the 
Froude number Fr = uj \/gh cos is small near the crit- 
ical angle or angle of repose for both the sand and the 
glass bead flows, it is natural to take the limit Fr — > 
in Eqs. (0E1), which suppresses the LHS of Eq. ©. We 
thereby obtain with some work an approximate equation 
for h: 



dh 
~dt 



■ a(h) 



dh 



d 2 h 



a(h) = yj gh 
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'dx " dx 2 ' ~ v " / v \2^h s (9) 2 

(4) 

where the "viscosity" v can be computed as v ~ dy/gh 
[2^ . This equation has solutions similar to those of 
Burger's equation, with a Burger's shock smoothed by 
the influence of the viscosity term. Thus, there is a so- 
lution consisting of a single hump propagating down the 
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slope with velocity a(/i), with a smooth structure deter- 
mined by the competition between this nonlinear velocity 
term on the LHS of Eq. (@J and the viscosity term. 

Turning to the full system of Eqs. (|2l3j) , however, we 
observe a potential flaw in this approach [l£j. This full 
system is hyperbolic with characteristic velocities 



a ± 




(5) 



If the velocity appearing in Eq. (HJ) does not obey a < c+, 
then Eq. (0} predicts a structure that moves faster than 
the maximum rate at which information can be propa- 
gated in the full system of equations, which is clearly 
impossible. In these circumstances, the Burger's type 
solution transforms itself into a truly discontinuous solu- 
tion traveling at velocity c + which is described by 
the full system Eqs. (2-3) rather than by Eq. gj). 

Using the rheology determined by Eqs. (QJ, one can 
see that for SGB, taking a value of a = 4/3, this dis- 
continuous solution will only develop for heights above 
h/h s ~ 6, whereas for RNSG, it will develop as soon 
as the avalanche has a height h/h s « 1.3 (the precise 
threshold depending on the value of 6). This quanti- 
tative height condition is consistent with the data in 
Fig. 03 where h m /h 3 < 6 for SGB and h m /h 8 > 1.3 
for RNSG. Further, calculation of the experimental ratio 
a/c + using the rheology of Eqs. (QJ yields a/c + < 0.75 
for SGB and a/c + > 1.1 for RNSG, which is also consis- 
tent with this picture. Thus, we conclude that the glass 
bead avalanches reflect smooth solutions of Eq. (J3J, with 
a < c + , whereas the sand avalanches represent discontin- 
uous solutions of the full system, traveling at velocity c+. 
The latter avalanches propagate into a quiescent bed be- 
cause they are traveling at the characteristic velocity for 
the medium. The glass bead avalanches are analogous to 
"flood waves" in river flows, whereas the sand avalanches 
are analogous to "roll waves" in these flows [HI [23. Note 
that ahead of the flowing avalanche, the moving material 
propagates into a material at rest, which is presum ably 
in a state close to the critical Mohr- Coulomb state |28j . 
Unlike the flowing state, for which <j xx « a zz , in this 
critical state <j xx > a zz . Thus the transition region, in 
which the flow accelerates from rest into a pseudo-steady 
state described by the continuum theory, can be viewed 
as a region of passive Rankine failure, through which the 
compressive stress parallel to the bed, a xx , is decreasing 
with time. The mechanics of this region is complex, and 
cannot be described by the Saint- Venant equations alone. 

Finally, we point out that in the linear theory of the 
instability of steady flows, developed by Forterre and 
Pouliquen, the criterion a < c + corresponds to the stable 
regime of these flows with respect to wave disturbances 
Md. Thus our observation of discontinuous avalanches 
for sand and smooth avalanches for glass beads dovetails 



nicely with their observation that steady flows of sand are 
far more unstable to such disturbances than are steady 
flows of glass beads. 
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